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ABSTRACT It is shown that there is a second properly nor- 
mahzed KNO scahng function, nPn{n/n) = f{z), which has 
certain advantages in the analysis of KNO scaling. First, the 
nPn are not influenced by the statistical and systematic uncer- 
tainties of n hence f>{z) provides more selective power than the 
original KNO scaling function nPn{n/n) = t/^^z). Second, the 
new scaling function generates scale parameter a = 1 since it 
depends only on the combination of z and the scale parameter 
of i){z). An analysis of f{z) is given using e~^e~ annihilation 
data for charged particle multiplicity distributions. 
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One of the most influential contributions to the analysis of multiplicity distribu- 
tions was made more than 20 years ago by Koba, Nielsen and Olesen [1]. They 
put forward the hypothesis that at very high energies the probability distributions 
Pn for detecting n final state particles exhibit a scaling law of the form 

Pn = li-('i). (1) 



n \n 

That is to say, the nPn measured at different energies (i.e. n) scale to the universal 
curve t/j when plotted against the multiplicity n rescaled by the average multiplicity 
n. This is the famous Koba-Nielsen-Olesen (KNO) scaling hypothesis. Defining the 
scaling variable z = n/n the scaling function i(j{z) must satisfy two normalization 
conditions: 

fOO 

^{z)dz=l (2) 



/ 

^0 



and 

ziIj{z) dz = 1, (3) 



/ 

^0 



i.e. z = 1. Obviously, the moments Aiq oi'ilj{z) are independent of collision energy 
if Eq. (1) is satisfied. 



/•oo 

Mq= z'ii;{z)dz 
Jo 



_ (4) 
= zi = ni/n'^ = const. 

Concerning the KNO scaling form of Pn given by Eq. (1) let us make a fur- 
ther remark. For a mathematically specified distribution F{x) a is called a scale 
parameter if F{x) has the form [2] 

a>0. (5) 

The parameter a characterizes the dispersion, that how widely the values of x are 
spread around the typical value, reflecting thus the scale of the distribution. From 
Eq. (1) it is seen that in the KNO scaling form of Pn the average multiplicity 
serves as a scale parameter. This coincides with the fact that if KNO scaling holds 
the dispersion of P„ is measured, up to a constant of proportionality, by n. The 
form Eq. (5) is obeyed by the scaling function too because the requirement z = 1 
usually generates a scale parameter cr 7^ 1 for ip{z). 

In the phenomenological analysis of multiplicity distributions one of the suc- 
cessfully applied probability laws is the negative binomial [3] . It has the form 
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where n and k are free parameters and 

Vik,n) = kik+l)...ik + n-l) = ^^ (7) 

denotes the Pochhammer polynomial for (possibly non-integer) rising factorial 
powers [4]. Eq. (6) has a nice KNO scaling form: in the limit n/k ^ 1 and n/n 
fixed 

nPn -^ai^) (8) 

where 

t{;^{z)=Af-d''z''-^exp{-9z) (9) 

is the gamma distribution with normalization constant J\f = T~^(k), scale param- 
eter 9 and shape parameter k. if^aiz) may have a variety of shapes, depending on 
the value of k, whereas 9 reflects its scale, the degree of scatter of z around z. The 
moments Aig of the gamma distribution are given by 

Mg = 9-'^V{k,q). (10) 

From the requirement Aii = 1 expressed by the normalization condition Eq. (3) 
we see that the scale parameter 9 in Eq. (9) is enforced to satisfy 

9 = Mi\,^^=V{k,l), (11) 

i.e. 9 = k. The scale parameter coincides with the shape parameter of the KNO 
scaling function. KNO scaling holds for energy independent k. 

Since the work of Koba, Nielsen and Olesen testing the validity of the scaling 
hypothesis and the analysis of the scaling function '0(^) sixe of permanent interest 
in multiparticle phenomenology [5]. But there is a noteworthy fact that escaped 
attention: besides il^{z) there is a second properly normalized scaling function 
obeyed by the Pn- The requirement Eq. (3) for the first moment of i^{z) defines a 
KNO scaling function normalized to 1, (fi{z) = ziIj{z). This yields the scaling law 

nPn = ^(^^^ (12) 

for the multiplicity distributions, the nPn should lie on the universal curve f{z) 
if Eq. (1) holds. The moments of the two KNO scaling functions are related by 

Investigating nPn instead of nPn has the obvious advantage that the statistical 
and systematic errors of n do not give contribution to the experimental uncertainty 
in the shape of the scaling function. Therefore (p{z) can be more selective between 
various theoretical predictions than '0(^). Another noticeable property of ip{z) 
arises in connection with Eq. (5). Since the KNO function t(j{z) obeys this form, 
with a constrained by Eq. (3), the scaling law in Eq. (12) can be written as 

nPn = ^(^\ (13) 
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The nPn depend only on the scaling variable z rescaled by the (energy indepen- 
dent) scale parameter a of iIj{z). Thus (p{z/a) provides a properly normalized 
KNO scaling function with scale parameter 1 and unconstrained M.i. The gamma 
distribution exhibits the scaling law Eq. (13) with a = d~^. Using the new scaling 
variable 

w^ze= (14) 

— 1 

{6 being restricted according to Eq. (11)) one arrives at the scaling function 

(p{w) = M ■ exp(-M;) (15) 

with M = T~^{k). The nPn should fall onto the above scaling curve if the nPn 
follow a gamma distribution with constant shape parameter k. 

It is of interest to find the family of distributions obeying the form Eq. (5). On 
the basis of the results of Ferguson [6] let us consider the following three different 
types of probability laws on [0, oo) 

a) Lognormal distribution 

h) Pareto distribution, possibly inverse power distribution 
c) Generalized gamma or Stacy distribution. 

If X is a random variable distributed according to a), h) and c) above, the random 
variable y = \n. x will have a 

a') Normal distribution 

h') Exponential distribution, possibly negative or translated 
c') Generalized normal distribution 

respectively. According to Ferguson these distributions play an important role 
in statistics since they are the only known ones that generate scale parameter 
families Eq. (5) of densities admitting complete sufficient statistics when samples 
of size > 2 are available [6] . 

Among the six probability laws a) — d) we turn our attention to the gener- 
alized gamma or Stacy distribution. The ordinary gamma, Eq. (9), and thus the 
exponential law h') appear as a special case. Moreover a) — 6), the lognormal and 
Pareto laws (the former one being very important in multiparticle physics [7]) can 
be obtained from c) in certain limits of the parameters [6]. We summarize some 
basic properties of this distribution to be used later. Further details can be found 
in the textbook of Johnson and Kotz [8] . The generalized gamma distribution has 
the form 

{z)=U- e^z^"^-^ exp(-6'2'^) (16) 

with scale parameter 6 and shape parameter k as in Eq. (9) and normalization 
constant M — iJiV~^{k). The additional parameter, the exponent generalizes 
the gamma distribution in a simple and useful manner: Eq. (16) corresponds to 
a power with exponent l//x of a gamma variable. This can be seen most clearly 
from the structure of the moments: 



(17) 
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the qth moment of the generahzed gamma coincides with the (g//i)th fractional 
moment of the ordinary gamma at fixed 9 and k. To get properly normalized KNO 
function 'ip{z) having M.i = I the scale parameter 6 must satisfy 



If the nPn follow a generalized gamma distribution with energy independent shape 
parameter k and exponent /j, the nPn obey the scaling law expressed by Eq. (15) 
with the corresponding normalization constant J\f and scaling variable 



as can be checked through eqs. (17) and (18). This gives back Eq. (14) for ^ = 1. 
It is worth noticing that the exponent /x is absorbed entirely into the rescaled 
variable w. 

The above extension of the gamma distribution is not new in the physical lit- 
erature. For a large class of nonlinear stochastic processes with pure multiplicative 
noise Eq. (16) arises as the stationary solution of the corresponding Fokker-Planck 
equation [9] . On the basis of the previous work the generalized gamma distribution 
was reviewed by Carruthers and Shih [5] pointing out that numerous theoretical 
KNO scaling functions are special cases. In its restricted form with M.i = 1 and 
6 given by Eq. (18) V'ggI-^) rediscovered by Krasznovszky and Wagner [10]. 
In a series of papers they carried out a pioneering work in determining the mag- 
nitude of the exponent n for a large amount of data available in various reactions, 
see [11] and references therein. Here we consider briefly the arguments made by 
Dokshitzer about the higher-order perturbative QCD effects on the shape of 'ip{z) 
in e~^e~ annihilations [12]. 

Over the past years much attention has been focused on the description of 
multiplicity distributions for hard processes in the framework of quantum chromo- 
dynamics (for a recent review see ref. [13]). Although the calculations predicted 
KNO scaling, which is indeed the case in e+e" annihilations, the scaling function 
proved to fall off exponentially in contradiction to the experimental results. Ac- 
cording to observations [14] the tail of i^iz) decreases faster than exponential. In a 
recent attempt Dokshitzer was able to take into account more precisely the energy 
conservation in parton cascades by involving next-to-next-to-leading QCD effects 
in the determination of 'ip{z) and its moments [12]. The obtained KNO scaling 
function behaves according to 



e = Mi 



r^{k,i/fi). 



(18) 



w^z''e = 



II ( 2;/^+! - z^' ) 



(19) 



%l){z) oc ey.^{-[DzY) 



(20) 



for large z with calculable constant D and exponent 



1 



(21) 



11 = 



l-l{as{Q)y 
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In the above formula as is the QCD running couphng constant corresponding to 
the hardness scale Q and 

.K(a)) = ^ P2) 

is the multiplicity anomalous dimension responsible for the energy growth of the 
average multiplicity n. In ref. [12] it was concluded that the better account of 
conservation laws drastically reduces the higher-order correlations and thus the 
width of ip{z). The exponent fj, was found to be ft; 1.6 which agrees with the 
results obtained in [11]. 

The comparison of eqs. (20) and (16) suggests that the generalized gamma 
distribution is well suited to characterize how the higher-order corrections of per- 
turbative QCD affect the shape of 'il^{z). Let us take Eq. (9), the ordinary gamma 
with exponentially decreasing tail as a reference distribution. One can model the 
suppression of correlations by raising the gamma variable z to the power l/jJ, < 1. 
With increasing z this procedure reduces more and more the contribution of the 
exponentially falling tail to the multiplicity moments zi. The original distribution 
will be squeezed, decreasing faster than exponential and having fractional multi- 
plicity moments with rescaled rank q/ /j, in Eq. (10). The same method can be 
applied in the opposite case too: one can model possible enhancement of correla- 
tions by raising the gamma variable 2: to a power > 1. Recall that the exponent 
H in Eq. (16) has already been determined for numerous reaction types [11]. This 
can provide valuable information on the degree to which the correlations are en- 
hanced/suppressed relative to the KNO scaling form, Eq. (9), of the negative 
binomial distribution. 

Having derived a new KNO-type scaling law for the multiplicity distributions, 
Eq. (12), it is tempting to check its validity on experimental data. We analysed the 
charged particle multiplicity distributions in e~^e~ annihilations from the ALEPH, 
AMY, ARGUS, DELPHI, HRS, L3, OPAL and TASSO collaborations in the cm. 
energy range ^/s = 9.4 — 91 GeV [14]. Our aim was twofold. First, to determine 
how closely the nP^ from different energies follow a universal scaling curve f{z). 
Second, to measure how our decisive power changes, in terms of statistics, if 
(p{z) is fitted to theoretical predictions instead of il^{z). 

Figure 1 displays the experimental data for nPn plotted against n/n (since 
the points are densely populated no different symbols are used to represent differ- 
ent data sets). According to expectation, the nP^ follow a unique scaling curve 
in agreement with the observed KNO scaling behavior of the nPn [14]. Another 
expected feature in connection with (p{z) is the increase of fitting proce- 

dures. The nPn are not influenced by the statistical and systematic errors of n 
hence the shape of (f{z) is less uncertain than that of ij{z). Studying each data 
set separately we carried out fits both to nP^ and nP^ (by a theoretical f{z) 
and i(j{z) respectively) in order to measure the significance of this effect in differ- 
ent circumstances. We used Eq. (16) for the theoretical ip{z) with normalization 
constant 2J\f (because of the charged particle data), scale parameter 9 given by 
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Eq. (18) and exponent fixed at the value /j, = 1.6. According to our results 
the increase in is substantial: among the 11 analysed data sets it was found 
to be > 50 % for 9 distributions. In the fitted values of the shape parameter k 
no significant change was observed. The minimum and maximum increase in 
correspond to the ALEPH and TASSO 34.8 GeV data respectively. In the former 
case the variation is only marginal: 

X^d.o.f. = 11/23 and x^/d-o.f. = 13/23. 

In the latter case the deviation is much more significant: 

X^d.o.f. = 11/16 and x^/d-o.f. = 29/16. 

An approximately 60 % increase in x^ is obtained when the fitting procedure 
was performed on the 11 distributions simultaneously using variable exponent /U. 
Again, no significant change was observed in the values of the fit parameters. 
These are A; = 6.06 ± 0.34 and 11= 1.29 ± 0.03, the exponent /i measuring the de- 
gree of deviation from an exponentially falling tail proved to be smaller than in 
refs. [11,12]. The scaling function ip{z) corresponding to the fitted parameters is 
represented by the solid curve in Figure 1. 

Let us summarize our main results. We have demonstrated that besides nPn 
the more simple combination nPn also scales to a universal curve in the variable 
n/n if KNO scaling holds valid. This somewhat unexpected behavior follows from 
the fact that the second normalization condition for V'l^), Eq. (3), defines a second 
properly normalized KNO scaling function, f{z) = zil){z). Thus Eq. (3) provides, 
on the one hand, certain restrictions: the first moment of ip^z) is constrained to 
be All = 1 and this usually generates a scale parameter cr 7^ 1 for ipi^z). On the 
other hand, Eq. (3) defines a KNO scaling function without the above restrictions. 
The gth moment of ^{z) coincides with the {q + l)th moment of il}{z), further, 
^{z) generates scale parameter a = 1 since it depends only on the combination of 
z and the scale parameter of il^{z). Perhaps the most profitable feature of the new 
scaling function lies in the fact that the statistical and systematic uncertainties 
of n are not involved by the nPn- This makes ^{z) more selective than '4>{z). 
According to our experience gained in the analysis of 11 multiplicity distributions 
the increase in statistics if the nP^ are fitted to a model prediction instead of 
fiPn exceeds 50 % for the majority of distributions. On the basis of the above 
features we conclude that the analysis of the scaling function ip{z) may become 
a useful method in the exploration of physical processes that give rise to KNO 
scaling behavior for the multiplicity distributions. 
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Figure Captions 

Figure 1. The KNO scaling function nPn{n/n) = f{z) for charged particle mul- 
tiplicity data in e~^e~ annihilations. The marked points represent ALEPH, AMY, 
ARGUS, DELPHI, HRS, L3, OPAL and TASSO data in the cm. energy range 
y/s = 9.4 — 91 GeV. The solid curve represents the theoretical (p{z) corresponding 
to Eq. (16), see the text for details. 
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